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<d : Abstract 

. Mature spiral galaxies consist of stars, gases and plasma distributed approximately in a thin 

r~| I disk of circular shape. The rotation velocities quickly increase from the galactic center and then 

Ph' achieve a constant velocity from the core to the periphery. The basic dynamic behavior of mature 

O ■ spiral galaxies are well described by simple models balancing Newtonian gravitational forces against 

^ I the centrifugal forces associated with a rotating thin axisymmetric disk. 

^ ' The Freeman model assumes a mass density distribution decreasing exponentially with radius 

in a thin disk extending to infinity. This assumption of exponentially decreasing mass was inspired 
, by the observed exponentially decreasing luminosity curves. This exponential mass distribution 

I produces rotational velocity profiles which increase from the galactic center to a maximum toward 

the outer core (as observed), but then decreases out to the periphery. This assumed exponential 
fSj . distribution is unable to adequately describe the measured constant (flat) velocities from the outer 

C5 I core to the periphery. By contrast, Mestel's finite disk model assumes a mass distribution decreasing 

Tlj" ■ more slowly (inversely with radius) that results in a constant rotational velocity across the entire 

. disk. But this inverse radius dependence of the mass distribution has an unrealistic singularity at the 

I galactic center and does not describe the central core rotation properly. Thus we propose combining 

' the Freeman and Mestel disk models to utilize their strengths and eliminate their deficiencies. We 

utilize the Freeman formula for the central core region, and then the Mestel formula beyond the 
central core to the galactic rim to describe the mass density distribution in a thin axisymmetric 
disk. Employing a connecting scheme for the transition from the Freeman to Mestel formulas, 
this combined model produces rotation curves comparable to the measurements for mature spiral 
galaxies. 

For a more general thin-disk model, we develop an alternative computational method to solve 
for the mass density distributions for various given measured rotation curves, instead of calculating 
the rotation curves from assumed mass density formulas (as assumed by Freeman and Mestel). 
We employ a thin axisymmetric disk of finite radius, and then compute radial mass densities that 
balance the Newtonian gravitational forces against the centrifugal forces (due to galactic rotation) 
at each and every point. The computational solutions for typical rotation curves show corresponding 
mass density distributions which decrease approximately exponentially in the central core (similar 
to the Freeman formula), and then transition to a slower inverse radius decrease (similar to the 
Mestel formula) to the periphery. 

Thus these diverse approaches yield similar results and are mutaully self-consistent with New- 
tonian gravity/dynamics. 
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1 Introduction 

1.1 Observational Knowledge of Radial Galactic Rotation Profiles 

Telescopic images of mature spiral galaxies reveal most of the stars, gas and plasma reside in an 
approximately circular disk that is very thin compared with its radius (Refs.[T]-[2]). The data on 
galactic rotational velocity profiles (Refs.[3]-[8]) of mature spiral galaxies are characterized by a 
rapid increase from the galactic center, reaching a nearly constant velocity from the outer core to 
the outer periphery. This basic feature may be idealized as 

V{r) = 1 - e-'l^^ , (1) 

where V{r) denotes the dimensionless rotational velocity measured in units of a characteristic 
velocity Vq and r the radial coordinate from the galactic center. The parameter Rc is a description 
from the data of the various "core" radii of different galaxies. Typical galactic rotational profiles 
described by ([1]) are displayed in Fig 1. As indicated by the measurement data, the rotation velocity 
typically rises linearly from the galactic center (as if the local mass was in rigid body rotation), 
and then reach an approximately constant (flat) velocity out to the galactic periphery. 

The observed galactic rotation curves (Eq([T]) and Fig.l) can not be explained by simply applying 
the so-called orbital velocity law, derived for a spherically symmetric gravitational field which is 
applicable to the Keplerian rotation of our solar-planet system, but not disk-like galaxies. In fact, 
the galactic mass distribution calculated by the orbital velocity law applied to these constant (flat) 
galactic rotation curves yields an increasing mass with radius, contrary to the measured galactic 
luminosity curves which decrease exponentially with radius. 



1.2 Thin-Disk Gravitational Models 

For a thin rotating galactic disk (revealed by telescopic observations), a balance can be achieved 
between the Newtonian gravitational forces and centrifugal forces at each and every point. Because 
the gravitational field of a thin disk is not spherically symmetric, the orbital velocity law is not 
applicable. As presented by Feng & Gallo [9], the equation of force balance in an axisymmetric 
thin disk can be written as 



2n 



(r cos (j) — r)d(j) 
(f^ _j_ _ 2rr cos (j))'^/'^ 



V{r 



i2 



p{f)hfdf + A = 0, (2) 



where all the variables are made dimensionless by measuring lengths (ex., r, r, h) in units of the 
outermost galactic radius Rg, disk mass density {p) in units of Mg/R^ with Mg denoting the total 
galactic mass, and rotational velocity [^(r)] in units of the characteristic galactic rotational velocity 
Vq. The disk thickness h is assumed to be constant and small in comparison with the galactic radius 
Rg. Actually the physically meaningful quantity here is the combined variable (ph) that represents 
the effective surface mass density on the thin disk. As long as the disk thickness h is much smaller 
than Rg, its mathematical effect is inconsequential to the value of (ph). The gravitational forces 
of the finite series of concentric rings is described by the first term (double integral) while the 
centrifugal forces are described by the second term. 

In ([2]), we call the dimensionless paremeter A "galactic rotation parameter", given by 

where G denotes the gravitational constant, Rg is the outermost galactic radius, and Vq is the 
characteristic velocity (which is often chosen as the maximum asymptotic rotational velocity). 
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According to ([2]), the corresponding rotational velocity V{r) can be determined from a given 
distribution of p{r) and the value of A, or vice versa. 

Due the complexity of the equation, there are limited analytical formulas for a few simplified 
cases such as the Freeman disk and Mestel disk. Unfortunately neither of these simplified equations 
can satisfactorily describe the complete rotation curve given by ([1]). 

By contrast, the numerical method presented by Feng & Gallo [9] solves an equation system 
including ([2]) and a conservation constraint for constant total mass of the galaxy Mg, e.g., 

27r [ p{r)hfdf = 1 , (4) 

to determine p(r) and A from a given V{r). Actually the same computational code can also be 
conveniently used to determine V{r) and A from a given p{r), with a slight modification. Thus, 
the methodology of Feng & Gallo [9] represents a more general approach with the axisymmetric 
thin-disk model to describe the galactic rotation for any measured rotation curve. 

2 Thin-Disk Models for Predicting Rotation Curves from Given 
Mass Distributions 



To facilitate a unified treatment for various thin-disk models, we discretize the governing equations 
([2]) and @ by dividing the one-dimensional problem domain [0, 1] into a finite number of line seg- 
ments called (linear) elements. As described by Feng & Gallo [9], each element covers a subdomain 
confined by two end nodes, e.g., element i corresponds to the subdomain [rj, rj+i] where and rj+i 
are the nodal values of r at nodes i and respectively. With each of the N — 1 elements mapped 
onto a unit line segment [0,1] in the ^-domain (i.e., the computational domain), independent 
residual equations can be obtained from the collocation procedure, i.e., 

7V-1 



E 



=1 



E{mi) K{mi) 



dr 1 

p{Ohm-^-^di + -AV{nf = G, (5) 



where K{m) and E{m) denotes the complete elliptic integrals of the first kind and second kind, 
with 

"^^(^)"F^rz\- (6) 

The residual equations ([5]) can be used to compute either the A^ nodal values of V{ri) from 
given distribution of p{ri) or the distribution of p{ri) from a given set of V{ri), with given values of 
A and h. Without loss of generality, the value of h is assumed to be 0.01 as comparable with that 
observed for the Milky Way galaxy. If the constraint equation Q is also used with a discretized 
form 

df 



2^E 



pmm:77d^-^ = o, (7) 

the value of A can also be determined as part of the numerical solution. 

These generally applicable equations are conveniently used for computing variables, even when 
analytical formulas are available for some special cases. Hence, a unified treatment for all cases is 
established for convenient comparison and analysis. Moreover, as discussed by Feng & Gallo [9j, 
imposing a boundary condition at the galactic center r = for continuity of derivative of p, i.e., in 
discretized form 

pin) = p{r2) , (8) 
is desirable for obtaining high-quality numerical solutions consistent with reality. 
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2.1 Freeman's Disk with Assumed Exponentially Decreasing Mass Distribution 
with Radius 

Freeman's (Ref.|10j) exponential disk gravitational model assumes a surface density function given 
by 

p{r)h = e-'/^'^ , (9) 

where Rd is a dimensionless scale length measured in units of Rg. Freeman assumed the disk 
extends to infinity, so the solution can be written as 

AV{rf = AirRay^ [h{y)Ko{y) - h{y)Ki{y)] , (10) 

where y = r/{2Rd) and Iq, Ii, -fCo) -f^i are modified Bessel functions. 

Figure [2] shows the velocity V{r) calculated with the present computational code, based on ([5]), 
for an exponential disk truncated at r = 1 with R^ set at 0.2 and A at 1.0, and points marked as 
"x" calculated with (llOp . The comparison is generally quite good especially in the central core. 
The slight discrepancy around the galactic rim r = 1 is expected due to the difference between 
truncated and infinite disks. A slight difference between the cases with = 101 and 401 indicates 
that the numerical discretization effects on the computed results are insignificant here. 

If the velocity profiles V{r) shown in Figure [2] for N = 101 and 401 are used as known inputs, 
the disk mass density p[r) can be computed by solving 

fir f'^ 
"? Jo 

For Rd = 0.2, 

/ e-^^fdf = 0.038383 . 
Jo 

With a boundary condition, 

p{r2) - /)(0) = ra , 

imposed at r = 0, the computed disk mass density profiles, corresponding to the given velocity 
profiles V{r) with N = 101 and 401, both recover that given by ([9]) with Rd = 0.2. Moreover, the 
computed Mg = 0.241167 (ss 2 x vr x 0.038383), and A = 0.999768 and 0.999987, are respectively 
verified for N = 101 and 401. Thus, the correctness of the computational code used herewith is 
verified. 

The Freeman disk with an exponentially decreasing mass density cannot produce a constant 
rotational velocity outside the central core, although the increasing velocity in the central core is 
reasonably described. 




2.2 Mestel's Disk with Assumed Mass Distribution Decreasing Inversely with 
Radius 

Another analytical model is Mestel's (Ref.|llj) disk, which has an assumed mass density distribution 
that decreases inversely with radial position. The further assumption of finite disk yields the 
following equation with a termination function written as 



p{r) 



A 



2'Khr 



1 sin (r) 

vr 



:ii) 
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to yield a constant velocity profile V{r) = 1. However, the Mestel's disk has a mass density 
distribution becoming singular at the galaxy center r = 0. This mathematical singularity at the 
galactic center is circumvented by setting the nodal value of p{ri) equal to the value of p{t2), 
equivalent to applying the boundary condition ([8|). With the rest of the nodes all given the p- 
values according to (fTTI) for A = 1.572914 (so that ([7]) is also satisfied), almost all of the computed 
nodal values of V{r) based on ([5|) are numerically equal to unity except for a few nodes around 
r = with significantly reduced values due to the boundary condition ([8]). 

This singularity problem is one undesireable feature that has limited the usefulness of Mestel's 
disk. And Mestel's incorrect description of the measured rotational profile in the galactic core is 
another reason for the limited applicability of this model. 



2.3 A Model Based on Combination of Freeman and Mestel Disks 



To reasonably capture the basic characteristics of an idealized galactic rotation curve (Eq.([T])), we 
combine the Freeman and Mestel disks, with Freeman describing the central core (r < R^) and 
Mestel outside the core (r > R^. In doing so, we assume that 



p(r) 



where 



and 



Rd 



Poe 
A 
27rhr 



+ 



-r/Rd 



1 sin (r) 

vr 



< r < iic 
,Rc<r<\ 



(12) 



Rc vr^l-i?2[i_2sin-i(i?c)/vr] 



A 



2vr/i Rre-^'-/^d 



1 --sm-\Rc 
vr 



so that both p and dp/dr are continuous at r = R^. Moreover, the mass conservation constraint 
([7]) can also be satisfied by setting the galactic rotation parameter as 



A 



N~l „i 



n=l 



df 



2vr ^ I p{C)hm^^dC 



Figure [3] shows the rotation curves V{r) for the Combined Freeman/ Mestel disks at several 
values of Rc- For Rc < 0.6 we find that V{r) increases from r = and approaches an asymptotic 
constant value in agreement with most galactic data. But the velocity profile of Freeman disk 
alone dominates as Rc approaches 1 (e.g., Rc = 0.8). The mass density profiles corresponding to 
the rotation curves in Figure [3] are shown in Figure HI Freeman's exponential mass density curve 
dominates in the interval when Rc > 0.5, except near the galactic rim r = 1 where the Mestel's 
formula takes over to force mass density to zero. 

Pronko (Ref. jl2j) arrived at a model similar to the Freeman/Mestel combination from a different 
perspective. 

To summarize, the Freeman/Mestel combination disk yields reasonable descriptions of galactic 
rotation velocity profiles. However, with assumed mass distributions even the combination lacks 
the flexibiility and versatility of the Feng/Gallo model (Ref.|i)J) (Eqs.2-4) which employs a compu- 
tational search technique to determine the mass distribution that satisfies any measured rotational 
velocity profile accurately. 
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3 Comparison with Thin-Disk Model for Numerically Predicting 
Mass Distribution from Rotation Curve 

The rotation curves in Figure [3] are not exactly the same as that described by ([1]) for a given value 
of Rc simply because of the differences in mathematical definitions. However, most of the curves 
are comparable. For example, at Rc = 0.015 a mass distribution described by (jl2p can produce a 
rotation curve fairly close to that by ([T]) as shown in Figure O The corresponding mass density 
distributions are shown in Figure [6l As also found by Feng Sz Gallo [9] , the mass density tends to 
decrease steeply from the galactic center, but beyond Rc, the mass density distributions decrease 
rather slowly towards the galactic periphery. This is consistent with the observed exponentially 
diminishing luminosity data, since both the temperature and opacity /emissivity are lower towards 
the periphery and we would expect the light distribution to decrease more quickly than that of mass 
towards the periphery. We believe that the assumption (often invoked by others) that the mass 
distribution exactly follows the light intensity distribution is fundamentally inaccurate due to lack 
of sound physical basis. Our computational method for solving p{r) based on observed V{r) avoids 
that inaccurate assumption and yields rationally reasonable mass distributions corresponding to 
the measured galactic rotational profiles. 

Although the combined Freeman/Mestel disk model can produce reasonable rotation curves 
comparable to our idealized curves (Eq(T]) especially for small Rc (Figures [5] and [6|) , its capability 
is limited due to the presumed form of mass density distribution. In contrast, our Feng & Gallo 
[S] model for numerically solving for p{r) from a given V{r) is much more versatile. For exam- 
ple, measurements usually reveal a more complex wiggling rotation curve, as shown in Figure [71 
mathematically given by 



The mass density distribution described by (jl2p cannot produce this more complex velocity profile. 
But from the computational Feng & Gallo [9] model, the corresponding wiggling mass density 
distribution (Figure [8] and Eq. (ll3p ) is computed by numerically solving ([5]) and ([7]). The value of A 
for the wiggling mass density curve is computed as 1.5770, as comparable to the reference value of 
1.5703. When solving for p(r) from a given V{r), however, we need a much larger (e.g., N = 4801 
used in the present work) than that for computing V{r) from a given p(r) (e.g., N = 401). This 
larger number of nodal points is needed for obtaining a good quality (smooth) curve with these 
more sharply varying input curves. 

4 Total Galactic Mass 

The measured rotational velocity profiles V{r) includes knowledge of Vq and Rg. With the computed 
value of the galactic rotation parameter A, the total galactic mass Mg can be calculated as 



To check viability we investigate the idealized rotational velocity profile V{r) of our own Milky Way 
galaxy shown in Figure[7]with Rc = 0.015. From the measurement data, the parameters appropriate 
for the Milky Way galaxy are Rc = 0.015, Vq = 2.5 x W^{m/s), and Rg = lO^(hght-years) = 
9.46 X 10^0 (m). We obtain ^ = 1.57 for the mass density distributions shown in Figure [8l and then 
the total galactic mass of the Milky Way is calculated from (jl4p as 



y(r) = 0.85 l-e-"/°-' 



+ 0.2e-^'' sin(67rr) + 0.2Vr 



(13) 



Mg = 




(14) 



AG 



Mg = 5.65 X W^^{kg) = 2.8 x lO" (solar-mass) . 
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This value is in very good agreement with Milky Way star counts of 100 billion, and further 
considering there is additional dust, grains, lumps, gases and plasma in all galaxies. In a subsequent 
companion publication, we shall employ Bulge+Disk models which are more compatible with visual 
observations of spiral galaxies and yield higher values of total galactic mass Mg. However, we 
emphasize the essential physics of galactic rotation is gravitationally controlled by the ordinary 
baryonic matter within thin galactic disks. 

Exploring solutions for a wide range of Rc values, Feng & Gallo [9] noticed that the computed 
values of A are typically within a small range around 1.6. 

5 Ordinary Baryonic Matter versus Dark Matter 

The successful models (Freeman/Mestel Combo and Feng/Gallo) described in this paper explicitly 
assume that all galactic matter is contained within the galactic disk. And these models implic- 
itly assume this is ordinary baryonic matter within the disk. Furthermore, the Freeman/Mestel 
Combo and Feng/Gallo models reproduce the measured galactic radial velocity profiles. Many 
other publications (Refs.p!2] - [35]) employ similar computational techniques based on Newtonian 
gravity /dynamics (or General Relativity) similar to Feng/Gallo and also find radial mass distri- 
butions of ordinary baryonic matter within the galactic disk that satisfy the measured rotational 
velocity profiles. These mass distributions decrease roughly exponentially from the galactic center 
in the central core, but then decrease less rapidly (inversely with radius) towards the periphery. 
This decrease is slower than the measured light distribution. Thus there is ordinary baryonic mat- 
ter within the galactic disk distributed towards the cooler periphery with lower emissivity/opacity 
and therefore darker. There are no mysteries in this rational scenario based on verified physics. 

However, there are many other publications (Refs.^36]-|57j) that take a different track. Anal- 
ogous to the measured radial light distribtuion, they start with a radial mass distribution within 
the disk that decreases approximately exponentially with radius. This is called the Mass/Light 
ratio assumption. But these models do not yield the measured rotational velocity profiles. To 
compensate, speculations are invoked re the existence of "massive peripheral spherical halos of 
mysterious Dark Matter" around the galaxies. But no significant matter has been detected in this 
untenable unstable gravitational halo distribution. This speculated Dark Matter is "mysterious" 
since it does not interact with electromagnetic fields (light) nor ordinary matter except through 
gravity. This Dark Matter must have other abnormal (non-baryonic) properties to maintain its 
peripheral spherical shape against the galactic rotation and gravitational attraction of ordinary 
matter. Furthermore, the assumption of a simple proportional relationship between mass and light 
intensity is not based on convincing physical principles. It does not accurately describe the the 
strong infiuence of both temperature and opacity/emissivity on light intensity. These deficiencies 
become apparent from edge-on views of spiral galaxies with a dark galactic line against a bright 
galactic background, revealing the substantial radial temperature and opacity/emissivity gradients 
across the galaxy. The concept of mass/light ratio is too inaccurate and oversimplified to apply to 
a complex galactic structure. Many unverified "mysteries" are invoked as solutions to real physical 
phenomena. 

To explain the failures associated with the assumption of the mass/light ratio, a different spec- 
ulation has been proposed (Ref.[58]) re deviations from Newtonian gravtiy/dynamics. But there is 
no independent experimental evidence of such deviations. Our use of Newtonian gravity/dynamics 
with sound computational techniques has proven sufficient to explain the observed rotational ve- 
locity curves. 

In conclusion, our approach utilizing Newtonian gravity/dynamics to computationally solve for 
the ordinary baryonic mass distributions within the galactic disk simulates reality and agrees with 
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observed data. 

6 Limitations and Strengths of Thin Disk Model 

Our simple thin disk model does not address many important features such as spiral structure, 
plasma effects, galactic formation, galactic evolution, galactic jets, black holes, relativistic effects, 
galactic clusters, etc.. 

It is well known (Ref. [1] ) that the internal gravitational behavior of a thin disk is much different 
than a sphere. This distinctly different behavior enables our thin disk model to describe the 
rotational dynamics of mature spiral galaxies, and the total galactic mass. 

Our thin disk model has finite radial extent. Beyond the galactic radius, we assume the density 
has dropped to the inter-galactic level, which is approximately spherically symmetric and thus no 
longer affects the galactic dynamics. We mention this because some others (Ref.|l]) have taken the 
relevant integrals to infinity, which we think is inappropriate. 

In our thin disk approach, we balance the gravitational forces against the centrifugal forces at 
each and every point. Thus, our solutions for the mass distributions and total galactic mass satisfy 
the measurements and ensure stability within the same context as similar calculations for our Solar 
System and Earth satellites. Some previous authors obtain solutions that are not gravitationally 
stable because they obtain incorrect mass distributions and incorrect galactic masses and do not 
satisfy the measured rotational profiles. Thus, their solutions are unstable, whereas our solutions 
are stable within the Newtonian context. 

Plasma effects are certainly active in the formation and evolution of galaxies from the original 
hot plasma (Refs.[59]-[61]). However, for the mature spiral galaxies we are addressing, the free 
plasma density has dropped to levels sufficiently low that plasma does not affect the predominantly 
gravitational galactic dynamics. This is evidenced in our own Solar System in which gravitational 
dynamics dominate even in the presence of solar wind, coronal mass ejections, comet tails, etc. The 
plasma in our Sun is stabilized by gravitational forces, even though plasma effects are very active 
within the Sun itself. Since our Solar System is approximately 1/3 distance from our Milky Way 
galactic center, we have our Solar System evidence for the dominance of gravitational forces within 
our own Milky Way galaxy, at least out to these radial distances. 

Again, we repeat, our thin disk model is sufficient to describe the rotational dynamics of mature 
spiral galaxies and the total galactic mass. 

7 Conclusions 

The measured rotation velocity profiles of mature spiral galaxies are successfully described with 
an appropriate combination of Freeman and Mestel thin disk gravitational models. These models 
start with assumed analytical radial mass distributions. 

By contrast, our approach utilizes Newtonian gravity/mechanics to computationally solve for 
radial mass distributions that satisfy the measured rotational velocity profiles. 

All these results reproduce the measured rotational velocity profiles, even though the funda- 
mental approaches are quite different. 

Also, all the radial mass distributions within the galactic disk are also mutually compatible. We 
compute that in the central galactic core, the mass density decreases approximately exponentially 
(similar to Freeman model). But beyond the core out to the periphery, we compute that the mass 
density decreases more slowly (inversely with radius similar to Mestel model). It is important to 
realize this mass distribution toward the galactic periphery is darker due to the lower temperature 
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and lower opactiy/emissivity in these outer regions. This is apparent from edge-on views of galaxies 
which display a dark disk-line against a much brighter galactic halo. 

Most other previous research assumes a galactic density decreasing exponentially with radius 
out to the galactic periphery, analogous to the measured light distribution. But this assumption 
(by others) is inaccurate since both the temperature and opacity /emissity are important but ig- 
nored variables. There is no simple relationship between mass and light. These prior models do 
NOT describe the measured velocity profiles, and speculations are invoked re halos of mysterious 
Dark Matter or gravitational deviations to compensate. The Dark Matter must have "mysterious" 
(non-baryonic) properties because there is no evidence of its existence and it is not responding to 
gravitational, centrifugal and electromagnetic forces in any known manner. By contrast, our results 
indicate no massive peripheral spherical halos of mysterious Dark Matter and no deviations from 
simple gravity. Our total galactic mass determinations are also in agreement with data. 

The controversy is summarized as follows. 

We believe there is ordinary baryonic matter within the galactic disc distributed more towards 
the galactic periphery which is cooler with lower opacity /emissivity (and therefore darker). 

Others believe there are massive peripheral spherical halos of mysterious Dark Matter surround- 
ing the galaxies. 
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Figure 1: Galactic rotational velocity profiles V{r) idealized from measurements according to ([T]) 
for Rc = 0.015, 0.03, 0.05 and 0.1. 
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Figure 2: Rotational velocity V{r) computed with ([9]) for a given surface density distribution 
p{r)h = e~^^ with N = 101 and 401 are represented by the solid curves. The data points (x) are 
calculated from Freeman's analytical solution (jlOp for an infinite disk with exponentially decreasing 
mass distribution. The computational codes are verified. 
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Figure 3: Velocity profiles V{r) corresponding to to the Combined Freeman/Mestel disks mass 
density distributions at different values of Rc witfi N = 401. The values of A for Rc = 0.01, 0.1, 
0.5, and 0.8 are computed as 1.5775, 1.6424, 1.9304, and 1.5096, respectively. 
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Figure 4: Selected density profiles p{r) of the Combined Preeman/Mestel disks corresponding to 
the velocity profiles shown in Figure [3] at different values of Rc with N = 401. 
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Figure 5: Reasonable comparison of V{r) for Rc = 0.015 corresponding to p{r) described by 
Combined Freeman/Mestel Disks (|12p and that computed from Feng &: Gallo [9J model (EqU]). 
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Figure 6: Reasonable comparison of p{r) for Rc = 0.015 for Combined Freeman/Mestel disks (jl2p 
(dotted curve) and that computationally predicted for the given rotation curve described by ([T]) 
(solid curve) and Feng & Gallo ^ model. 

The mass density decreases steeply from the galactic center, but then more slowly toward the 
periphery. 
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Figure 7: The more complex measured hypothetical wiggling rotation curve V{r) described by (jl3p 
together with that at Rc = 0.015 given by ([1]) (as labeled) as a smoother idealized comparison. 
Actually, these curves are similar to the measurements of our own Milky Way Galaxy, and we use 
this information to calculate the total galactic mass of the Milky Way Galaxy in the text. 
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Figure 8: The more complex wiggling mass density distribution from the Feng & Gallo [9] model 
corresponding to the hypothetical wiggling rotation curve V{r) (jl3p compared with the smoother 
rotation curve given by ^ for Rc = 0.015. 
The versatility of the Feng & Gallo [9] model is apparent. 



